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Abstract. We show that the authors of the commented paper draw their conclusions 
from the eigenvalues of truncated Hamiltonian matrices that do not converge as the 
matrix dimension increases. In one of the studied examples the authors missed the 
real positive eigenvalues that already converge towards the exact eigenvalues of the 
non-Hermitian operator and focused their attention on the complex ones that do not. 
We also show that the authors misread Bender's argument about the eigenvalues of 
the harmonic oscillator with boundary conditions in the complex- a; plane(Rep. Prog. 
Phys. 70 (2007) 947. 



1. Introduction 

In a recent paper Bowen et al [I] discussed the spectra of a class of non-Hermitian 
hamiltonians having PT symmetry. They calculated the eigenvalues of truncated 
matrices for the non-Hermitian Hamiltonian operators in the basis set of the 
eigenfunctions of the harmonic oscillator and argued that their results did not agree 
with those of Bender and Boettcher [2]. They concluded that the discrepancy may be 
due to the fact that the WKB method used by the latter authors is unsuitable for such 
problems. For example, they stated that "It is certainly not obvious that the physical 
conditions of the Bohr-Sommerfeld procedure should be valid for this non-physical path. 
The significant differences in the spectrum studied in this paper suggests that it is not 
valid" and also that "It is not clear whether the motion along paths in the complex plane 
has any physical significance for quantization". Curiously, the autors did not appear 
to pay attention to other methods for the calculation of the eigenvalues of those PT- 
symmetric Hamiltonians. For example, the WKB results were confirmed by numerical 
integration based on the Runge-Kutta algorithm [2|]3] as well as by diagonalization of a 
truncated Hamiltonian matrix in the basis set of harmonic-oscillator eigenfunctions [2] 
(a more detailed description of this approach was given in an earlier version of the 
published paper [6j). In addition to it, Handy [HE] obtained accurate upper and lower 
bounds from the moments equations. 
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The results, conclusions and criticisms of Bowen et al pQ are at variance with all 
what has been established after several years of study in the field of non- Hermit ian 
PT-symmetric Haniltonians [3]. The purpose of this comment is to analyse their 
calculations to verify if such criticisms are well founded. In section |2] we briefly review 
the diagonalization method used by the authors, in section [3] we analyse some of the 
models used by the authors to draw their conclusions; finally, in section H] we summarize 
the main results and draw our own conclusions. 

2. The method 

Bowen et al [1] calculated the eigenvalues of the class of Hamiltonian operators 



for s = 1,-1, i and N = 2,3,4,6,8. They resorted to matrix representations of the 
operators in the basis set of eigenfunctions {\n) , n — 0, 1, . . .} of the harmonic oscillator 
(s = 1, N = 2) and diagonalized truncated Hamiltonian matrices H^ M ^= (H mn ) m ~^ , 
where H mn = (m\H\n), for each of those cases. If the eigenvalues of the truncated 
matrices converge as M increases then the limits of those sequences approach the 
eigenvalues of the operator (JI|) . 

The characteristic polynomial for the matrix H^ M ^ will exhibit M roots W^ M \ 
n = 0, 1, . . . , M — 1. In the case of Hermitian operators (for example, s = 1 and iV 
even) all those roots are real because the matrix is Hermitian. It is not the case for the 
non-Hermitian operators. Furthermore, in the case of the Hermitian operators we know 
that the eigenvalues of the matrix approach the eigenvalues E n of the operator from 
above W { n M) > W^ M+l) > E n . On the other hand, there is no such variational principle 
in the case of non-Hermitian operators. Obviously, one has to be very careful when 
applying the diagonalization method (DM from now on) to non-Hermitian operators. 
In what follows we discuss some examples. 



H = p 2 + sx 



(1) 
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3.1. Case N = 2 

When s = 1 the matrix H is diagonal and yields the eigenvalues of the harmonic 
oscillator exactly: E n = 2n + 1, n = 0,1,.... On the other hand, for s — — 1 
the eigenvalues of the matrix do not converge as M increases. Therefore, they are 
meaningless and bear no relation with the eigenvalues of the Hamiltonian H = p 2 — x 2 . 
Surprisingly, the authors argued that "This Hamiltonian has a spectrum with odd 
symmetry about zero energy" . If one solves the eigenvalue equation with the appropriate 
boundary conditions in the complex- a; plane one obtains purely imaginary eigenvalues: 
E n = ±(2n + 

The authors went even further and stated that "Bender has also asserted that the 
spectrum of the simple harmonic oscillator (SHO) (r = 0) with a negative force constant 
has a discrete negative spectrum that is the negative of the positive force constant SHO; 
that is, E n = — hu(n +1/2)" . However, Bender [3] never draw such wrong conclusion. He 
discussed the harmonic oscillator H = p 2 + u 2 x 2 with eigenvalues E n (u) = (2n+l)u and 
eigenfunctions ip n {^>, x). The eigenfunctions behave asymptotically as ip n (u, x) ~ e - wx ' 2 / 2 
when \x\ — > oo. If we substitute — u for u the eigenvalues change sign but the 
eigenfunctions are no longer square integrable. However, if we rotate the variable 
7r/2 counterclockwise then the resulting eigenfunction ip n (cu,q) = ipn^—uj^iq) is square 
integrable. It is quite obvious that the substitution of —u for u changes the sign of the 
eigenvalues but the force constant (oc u 2 ) does not change. In fact, Bender [3] states 
that "Notice that under the rotation that replaces u by —u the Hamiltonian remains 
invariant, and yet the signs of the eigenvalues are reversed!" Therefore, it seems that 
Bowen et al [1] misread Bender's argument. 
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3.2. Case N = 3 

When s — 1 the eigenvalues of the truncated matrices do not converge as M increases. 
However, the authors state that "Here the spectrum was almost symmetric about zero..." 
in spite of the fact that the roots of the secular determinants are not valid approximations 
to the eigenvalues of the differential operator. 

The only interesting case is undoubtedly the PT-symmetric Hamiltonian operator 
for s = i. According to the authors "The calculation of the spectrum for the potential 
V = ix 3 yielded a complex spectrum". In this case the wedges in the complex- x plane 
where ip(x) vanishes exponentially as \x\ — > oo contain the real axis [2]. Therefore, 
one expects the DM to yield meaningful results. Our calculation shows that the 
complex eigenvalues of the truncated matrices do not converge as M increases, but 
there are real ones that certainly converge towards the results obtained by Bender and 
Boettcher [2] by means of the WKB method and numerical integration. In fact, Bender 
and Boettcher J2j[6] discussed the calculation of the eigenvalues by means of the DM. 
They concluded that the method is only useful when 1 < N < 4 and that the convergence 
to the exact eigenvalues is slow and not monotone because the Hamiltonians are not 
Hermitian. Table [1] shows the convergence of the lowest eigenvalues of the truncated 
matrices towards those obtained by means of the Runge-Kutta method (RK) and the 
WKB approach; they are real and positive as argued by Bender and Boettcher [6]. 

An interesting feature of the DM for non-Hermitian operators is that the 
characteristic polynomial of degree M does not exhibit M real roots as in the case 
of the Hermitian matrices. In the present case the truncated matrices also exhibit many 
complex eigenvalues but they do not converge as M increases. Another interesting 
feature is the behaviour of the approximate eigenvalues with respect to a scaling factor. 
Instead of using the eigenf unctions ip n ( x ) — (x\ n) of the harmonic oscillator p 2 + x 2 
we can try an alternative calculation with the scaled eigenfunctions a 1//2 ?/> n (ct£), where 
a is an adjustable scaling factor. In the case of Hermitian operators W^ M \a) exhibits 
a minimum because of the variational principle. On the other hand, in the case of the 
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complex potential V = ix 3 the approximate eigenvalue oscillates and exhibits a kind 
of plateau with oscillations of smaller amplitude. The optimal value of a is somewhere 
in this region. For example, we find that a ~ 1.4 is more convenient than the scaling 
parameter a = 1 used in the calculation shown in Table [U However, it is our purpose 
to show here only the results for the same basis set chosen by Bowen et al [lj. 

3.3. Case N > 4 

In these cases the DM only yields meaningful results for s — 1 and N even; that is 
to say: for the trivial Hermitian operators. When s = -1 or s = i the eigenvalues of 
the truncated Hamiltonian matrices do not converge and, consequently, they are not 
eigenvalues of the Hamiltonian operator. The failure of the naive application of the DM 
to the PT-symmetric cases is not surprising because the wedges in the complex-x plane 
where ip(x) vanishes exponentially as \x\ — > oo do not contain the real axis [2]. 

The PT-symmetric Hamiltonian operator H = p 2 — x 4 deserves special attention 
because several authors have proved that it is isospectral to the Hermitian one H = 
p 2 + 4a; 4 — 2x [3] (and references therein). Apparently, Bowen et al pQ were not aware 
of this relationship that could have convinced them that the former Hamiltonian does 
already have a positive spectrum. 

4. Conclusion 

It is clear that the discrepancy between the results of Bowen et al [1] and Bender and 
Boettcher [2] is merely due to the fact that the DM used by the former authors does not 
apply to some of the problems studied. Their conclusions were based on eigenvalues of 
the Hamiltonian matrices that do not converge. They only obtained meaningful results 
for the trivial cases of Hermitian Hamiltonians given by s = 1 and N even. In the 
only other selected case where the DM is expected to yield reasonable results, namely 
V(x) = ix 3 , the authors failed to find the converging real positive roots and simply 
focussed on the complex ones that do not converge. 
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Table 1. First eigenvalues of the truncated matrices of dimension M for V(x) 



M Eq E\ E2 -E-3 



10 


1.156101684 


3.73083496 








15 


1.156038818 


4.14942907 








20 


1.156383056 


4.109441589 








25 


1.156258544 


4.109537412 


7.553497517 




30 


1.156267013 


4.109170441 


7.562399797 


11.24884001 


35 


1.156266986 


4.10922^ 


3991 


7.562011977 


11.31452225 


40 


1.156267082 


4.10922^ 


3365 


7.562284307 


11.31372188 


45 


1.156267072 


4.10922^ 


3831 


7.562273020 


11.31452360 


50 


1.156267072 


4.10922^ 


3753 


7.562274330 


11.31442188 


55 




4.10922^ 


3754 


7.562273^ 


354 


11.31442413 


fin 




4.10922^ 


3753 


7.562273^ 


360 


11 01 AAO~\ IP. 


65 




4.10922^ 


3753 


7.562273^ 


354 


11.31442184 


70 




4.10922^ 


3753 


7.562273^ 


355 


11.31442182 


75 








7.562273^ 


355 


11.31442182 


80 












11.31442182 


RK 


1.156267072 


4.10922^ 


3752 


7.562273^ 


354 


11.314421818 


WKB 


1.0943 


4.0895 




7.5489 




11.3043 
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